On the Lifshitz tail in the density of states of a superconductor with magnetic 

impurities 



A.V. Balatsky* and S. A. Trugman 
Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 

(February 1, 2008) 

We argue that any superconductor with magnetic impurities is gapless due to a Lifshitz tail in 
the density of states extending to zero energy. At low energy the density of states v{E — > 0) 
remains finite. We show that fluctuations in the impurity distribution produce regions of suppressed 
superconductivity, which are responsible for the low energy density of states. 
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The role of impurities in superconductors is a rich sub- 
ject, going back to the pioneering papers by Abrikosov 
and Gor'kov jl| and by Anderson B. However the ma- 
jority of work has been concentrated so far on the "mean 
field" treatment of the impurity problem in supercon- 
ductors. Here we will address the role of the fluctuations 
of the distribution of magnetic impurities in an s-wave 
superconductor. 

It has been experimentally known for some time that 
the density of states (DOS) in a superconductor with 
magnetic impurities is far greater at low energies than one 
would expect from Abrikosov-Gor'kov theory |3). Using 
the suppression of the critical temperature to infer the 
pairbreaking parameter, one typically arrives at a sub- 
stantially lower DOS at E <C Ao than is observed. (Ao 
is the superconducting gap in the spectrum.) We suggest 
here that the observed deviations from the Abrikosov- 
Gor'kov theory at small E are caused by fluctuations in 
the impurity distribution and the Lifshitz tail in the DOS 
of an impure superconductor. 

We observe that for any, no matter how small, con- 
centration n of magnetic impurities in a superconductor 
there are fluctuations in the distribution of impurities 
across the sample. There are finite regions of high im- 
purity concentration, where the superconducting state is 
suppressed due to scattering. These large regions of es- 
sentially normal metal produce low lying, E <C Ao single 
particle states in the averaged density of states of the su- 
perconductor. It is clear that any singularity in the DOS, 
if one occurs, should be at E = due to the particle-hole 
symmetry of the superconducting state, which we assume 
here and which is preserved even with magnetic impuri- 
ties. We find that at low energy £ < Ao the DOS is: 

u(E) oc (l/A Lo )exp(-const L d ), (1) 

where v{E) is scaled with the normal state DOS, d is the 
dimensionality of space, Al is the mean level spacing in 
the fluctuation region of the size of length Lq = (CoO 1 ^ 2 ; 
where £o = nvf/Ao is the T = superconducting coher- 
ence length, and I is the mean free path. The constant 
in the exponent will be given below. The tail in the DOS 



of a superconductor is similar to the tail in the DOS of 
a semiconductor, the so-called Lifshitz tail M. 

For any particular model of impurity scattering (e.g. 
Born versus unitary scattering), we assume that there ex- 
ists a critical concentration n c at which a thermodynamic 
superconducting sample will become normal due to the 
pairbreaking effect of impurities. The specific value of 
n c obviously depends on the model. For the case of the 
Born scattering limit of magnetic impurities, within the 
Abrikosov-Gor'kov theory, n c = O(l)A N /(J 2 N$S(S+ 
I)), where Nq is the normal metal DOS, J is the mag- 
netic exchange between conduction electrons and the im- 
purity, and S is the magnitude of the impurity spin [[l]. 
This specific value is not important for our subsequent 
considerations. We will use n c as a model-dependent in- 
put to our final answer. All concentrations are given in 
terms of the dimensionless concentration per unit cell of 
linear size a. 
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FIG. 1. Fluctuation region of size L, with concentration of 
impurities n c inside the superconductor, is shown schemati- 
cally. The equilibrium concentration is n < n c . The fluctua- 
tion region has a metallic spectrum. Andreev reflection mod- 
ifies the spectrum of quasiparticles [^|. In any local probe 
of the DOS, e.g. an STM, one would find that the I - V 
characteristics have a gap in the outer region, but are gapless 
if measured at any point inside the fluctuation region. The 
average DOS of a superconductor as E — > hence will be the 
sample average of the DOS of the fluctuation regions. 

Here we will consider the case of arbitrary impurity 
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exchange strength. It is known that magnetic impurities 
induce intra-gap states ||. The energy of these states 
for large S is approximately ujq = A (l — J 2 NqS(S + 
1))/(1 + J 2 N$S(S + 1)). These impurity states have a 
wavefunction 4"(r) ~ exp(—r/^ Ua ) of size = £n(l — 
(wo/A ) 2 ) -1 / 2 > £o> where £n is the zero temperature 
superconducting coherence length. For subsequent con- 
sideration we assume that n c (^ /a) d » 1, generally true 
for realistic systems, so that intra-gap states are strongly 
overlapping in the region where the impurity concentra- 
tion is n c . Impurity states form an impurity band, cen- 
tered around loq ||. Fluctuations in the distribution of 
impurities lead to tails in this impurity band, which ex- 
tend to zero energy. 

Consider a fluctuation in the impurity distribution 
such that inside a region V(L) — L d 0, the local concen- 
tration of impurities is n c (averaged over distances much 
greater than £ Uo but smaller than L), as shown in Fig. 
1 . We assume that L 3> I > £ Wo , where I is the mean free 
path at the critical impurity concentration n c . 

The low energy single particle spectrum in the fluctu- 
ation region will be normal, since the local concentration 
is n c . The proximity coupling to the superconducting 
reservoir at the boundary cannot open up a gap at large 
distances ~ L 3> £o due to pairbreaking scattering. We 
ignore the region of size £n from the boundary where the 
gap is decaying. The single particle spectrum inside V(L) 
will be equivalent to the spectrum of a normal metallic re- 
gion with magnetic impurities in tunneling contact with 
a bulk superconductor. 

To verify that the spectrum of the fluctuation region is 
indeed gapless we have numerically calculated the spec- 
trum of a random superconductor in the mean field ap- 
proximation. Specifically, we considered the ID BCS su- 
perconductor with the Hamiltonian 

H = -t ^ C i,a C j,a + A * c i,T c U + h.C, + 
<i,j>,a i 

o a ,f)Ci.{) (2) 

i£V{L),a,0 

where i labels the sites of ID chain, V(L) are the impu- 
rity sites, t is the nearest-neighbor electron hopping, Aj 
is the pairing amplitude on the site i, J is the exchange 
coupling between the conduction electron and impurity 
spin, and S, is a random classical Heisenberg impurity 
spin on the site i. The last term in Eq. (||) describes 
the impurity scattering effects of the fluctuation region, 
which we assume to be in the middle of the supercon- 
ducting region. 

We consider a superconducting system of 40 sites with 
impurity spins present at a high concentration x on 10 
of these sites. This approximation was chosen to mimic 
the high impurity density fluctuation region, which is re- 
sponsible for the low energy DOS. For classical spins, the 
coupling J and impurity spin magnitude enter into the 



answer in the combination JS, hence the specific values 
of each of them separately does not matter. We have 
calculated the spectrum of quasiparticles in the mean 
field approximation, ignoring the self-consistency condi- 
tion for the gap S. The DOS for this model is shown in 
Fig. (2). 
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FIG. 2. The density of states is plotted for a ID BCS su- 
perconductor with 40 sites. There is an impurity region of 
10 sites, in which the superconducting gap A is taken to be 
zero. Classical Heisenberg impurity spins with random orien- 
tation occupy the impurity region with concentration x — 0.5. 
The solid line is for coupling JS = 0.1, and the dotted line is 
JS = 1.0. Other parameters are t = 1 and A = 0.5. The DOS 
is averaged over 1500 realizations. The fine-scale roughness 
in the middle of the band is due to finite size effects. 

Since A = in the impurity region, there are intra- 
gap states even for JS = 0. (There is only one such state 
for the parameters of Fig. 2.) We find that the intra- 
gap state evolves into an impurity band, and gradually 
fills the entire gap as the concentration or the coupling 
constant JS increases. This evolution of the impurity 
band is similar to the evolution of the band in doped 
semiconductors. The calculation confirms all the basic 
features one might expect: the appearance of impurity 
states inside the gap region, the growth of the impurity 
band, and finally the filling of states at low energies with 
nonzero ^(0). 
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FIG. 3. Same as figure (2), but with the superconducting 
gap A = 0.5 throughout the sample, including the impurity 
region. 

A similar calculation is shown in Fig. (3), but with 
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the mean-field superconducting gap A = 0.5 everywhere, 
including the impurity region. In this case, there are 
no intra-gap states for JS = 0. For small JS, intra-gap 
states first appear at ujq, which is just below the energy A 
of the uniform gap. As JS increases, the gap gradually 
fills in. A larger value of JS is required to completely 
close the gap than in Fig. (2), because the density must 
spread down from A rather than from the intra-gap levels 
that already exist in Fig. (2) at JS = 0. 

A similar problem for a metallic grain in the presense 
of time reversal violating fields (e.g. impurity spins) in 
contact with a superconductor was considered by Alt- 
land and Zirnbauer ||. At energies small compared to 
the Thouless energy Et = D/L 2 , one can ignore the spa- 
tially inhomogeneous solutions of the nonlinear-cr-model. 
In this limit the spectrum of the grain is given by ran- 
dom matrix theory. The single particle DOS in Ref. ||] 

is v L (E) = 1/A L (l + s " l 2^E%t L) ) ' wnicn S ocs to in- 
stant as E — > 0. Here A^ (not to be confused with the 
gap Ao) is the mean level spacing of the grain of linear 
size L, and vl{E) is averaged over all realizations of the 
random spectrum for grains of size L. We are interested 
in E <C cjo, where the constraint E <C Et is not im- 
portant since Et/ujq ~ (Ao/c^oXCoV-^ 2 ) an( i this ratio 
is small except in the limit ujq — » 0, where special care 
should be taken. We will not address this limit here. We 
believe that the result v(E) ~ const will still hold. 

If we make the assumption that the spectrum of the 
fluctuation region in Fig. 1 is equivalent to the spectrum 
of a normal metal grain, it is easy to estimate the average 
DOS v(E) from the distribution Pi(n c ;n) in the size of 
the fluctuation regions: 

y(E) ~ / ' dV{L)P L {n c -n)u L {E). (3) 

We now consider the probability distribution for a nor- 
mal region of volume V(L) with linear size L S> £o to oc- 
cur. This question is equivalent to finding the probability 
Pl(ji c ; n) of a fluctuation region of diameter L, taken to 
be spherical in d-dimcnsions, with a concentration of im- 
purities in this region equal to or greater than n c , while 
the average concentration is n. This probability can be 
easily evaluated, following, for example, the arguments 
of Refs. |||. We find: 

log(P L (n c ; n)) = 5<r = -V(L)<p(n c ; n) 

4>{n c \ n) ~ n c log(n c /n) - n c + n, (4) 

where 5a is the change in entropy due to a fluctuation 
with homogeneous concentration n c in the region V(L), 
and <j>{n c \ n ) is the entropy density for the discussed fluc- 
tuation, which is model dependent. Equation (g) ap- 
plies for small n and n c p| . Strictly speaking, Eq. (Q) 
gives the probability of a fluctuation with a concentra- 
tion equal to n c . In principle one should integrate this 



probability over the range n > n c to obtain the total 
probability that the normal region V(L) will occur. Tak- 
ing into account this effect will only change the coefficient 
in 4>(n c ; n) and the prefactor in Eq. (1). 

The ratio of the mean level spacing to the supercon- 
ducting gap is given by Al/A — n(n c , J, N )L~ d , where 
n{n c , J, No) is a model-dependent dimcnsionless function 
of n c ,J,N Q Q. With the aid of Eq. (g) and using A L 
we find: 

v{E)= dV{L)P L {n c ;n)v L (E) 

Jv(L ) 

~ Alhxp(-L d <p(n c ;n)) (5) 

This is our main result. In writing Eq. (ph the lower limit 
of the volume integration was taken at L = La = (tiol) 1 ^ 2 , 
when Et ~ Ao, because at smaller distances the gap 
acquires a nonzero mean value due to strong coupling to 
the bulk superconductor. 

The fact that v(E) is nonzero at arbitrarily small en- 
ergy implies that a superconductor with magnetic impu- 
rities is always gapless. This does not, however, mean 
that the system is not superconducting. A DC current 
can flow through the system (around the impurity re- 
gions) with no dissipation, i.e., there is a condensate. 
The dissipation is nonzero for essentially any AC current 
due to dissipation in the normal metal regions. 

A few comments are in order here, i) It should be 
noted that there is a qualitative difference between the 
DOS in the tails for a superconductor as compared to 
a semiconductor. In the case of a high impurity con- 
centration region of size L in a semiconductor, the en- 
ergy has a quadratic dispersion E — Eq = ir 2 /2mL 2 , 
where Eq is the lowest energy of the crystal composed 
of only the impurity atoms. This results in a v{E) oc 
exp(— const/ (E — Eo) d l 2 ) for the Lifshitz tail in a semi- 
conductor. The difference comes from the fact that in a 
semiconductor tails are formed near the bottom of the 
band, whereas in our case the destruction of the super- 
conducting gap leads to disordered normal regions, ii) 
The suppression of superconductivity occurs at quite a 
low concentration n c ~ 1%. This allows for substan- 
tial fluctuations of the impurity distribution inside V(L). 
However it is clear that the most important configura- 
tion responsible for the low lying states is the one with 
a nearly homogeneous distribution with local concentra- 
tion close to n c . Any fluctuations with local n(r) < n c 
arc ineffective for v{E — > 0). We expect any improve- 
ment of the above consideration will lead to corrections 
to (j>(n c ; n). iii) We have ignored the possible interactions 
between impurity spins. This does not have to be the case 
in real systems, where in order to suppress superconduc- 
tivity one has to have many impurities in regions of the 
size of the coherence length: n c (£o/a) d 3> 1. Interactions 
between spins in this situation may be important, as was 
pointed out by Larkin et. al. p3|| . iv) Similar consider- 
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ations should apply to any system with a spontaneously 
induced gap due to interactions, i.e. CDW, SDW systems 
and to unconventional, e.g. d-wave, superconductors. 

The present work is related to that of Larkin and 
Ovchinnikov jl4j]. They considered the DOS fluctuations 
for a disordered superconductor due to fluctuations of the 
gap Ao(r), and also found that the DOS is finite at small 
energies due to this process. We have considered here a 
different mechanism for generating a nonzero DOS. 
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Scalapino and M. Zirnbauer for useful discussions. We 
are grateful to G. Aeppli for bringing Refs.( §|) to our 
attention. Part of this work was done at ITP, Santa Bar- 
bara. This work has been supported by the Department 
of Energy and by NSF Grant PHY94-07194 at ITP. 
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